First-principles study of structural, elastic, and bonding properties of pyrochlores 
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Density Functional Theory calculations have been performed to obtain lattice parameters, elas- 
tic constants, and electronic properties of ideal pyrochlores with the composition A2B2O7 (where 
A=La,Y and B=Ti,Sn,Hf, Zr). Some thermal properties are also inferred from the elastic proper- 
ties. A decrease of the sound velocity (and thus, of the Debye temperature) with the atomic mass of 
the B ion is observed. Static and dynamical atomic charges are obtained to quantify the degree of 
covalency/ionicity. A large anomalous contribution to the dynamical charge is observed for Hf, Zr, 
and specially for Ti. It is attributed to the hybridization between occupied 2p states of oxygen and 
unoccupied d states of the B cation. The analysis based on Mulliken population and deformation 
charge integrated in the Voronoi polyhedra indicates that the ionicity of these pyrochlores increases 
in the order Sn-Ti-Hf-Zr. The charge deformation contour plots support this assignment. 

PACS numbers: 71.20.-b, 62.20.-x, 65.40.-b, 91.60.Ki 
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I. INTRODUCTION 

Materials with the pyrochlore (A2B2O7) lattice struc- 
ture have unique properties that make them ideal can- 
didates for applications ranging from high-permittivity 
dielectrics 1 , to ceramic thermal barrier coatings (TBC) 2 . 
potential solid electrolytes in solid-oxide fuel cells 3 -, or im- 
mobilization hosts of actinides in nuclear wasted. Many 
studies have tried to optimize the pyrochlore composition 
to obtain the most desirable properties for a particular 
application. One has to achieve maximum efficiency in 
a special property (high dielectric constant, low thermal 
conductivity, low migration energies, etc.), but at the 
same time satisfy other requirements such as physical 
stability, chemical or thermal compatibility with other 
materials involved in the design, etc. In many cases, no 
experimental information is available, or if there is infor- 
mation there is no agreement between different experi- 
ments, a consequence of the difficulties in growing pure 
samples. 

Therefore, simulations have become an ideal tool for 
systematic investigations of various properties as a func- 
tion of the chemical composition^iSiL2ii£. Most of the 
atomistic simulations studies are based in classical po- 
tential methods, where the interatomic interactions are 
parametrized with effective potentials fitted to reproduce 
experimentally known data, as the unit cell dimensions 
of a variety of oxide pyrochlores. The complexity of the 
chemical elements involved (mostly transition metals and 
rare earths) makes desirable the use of ab initio calcula- 
tions, but the large size of the unit cell in pyrochlores (88 
atoms) had restricted the application of this technique. 

The limitations of these classical models have recently 
been pointed out by Panero et. aliii in the study of 
the energetics of cation-antisite defects. The stability of 
pyrochlores and the resistance to amorphization by irra- 
diation has been correlated with the propensity of the 



ions to create these point defects^. With classical mod- 
els, the defect-formation energies show a dependence on 
the cationic radius: compounds with similar radii for A 
and B are more radiation resistant than compounds with 
very dissimilar radii, that show greater susceptibility to 
amorphization. The effect of the B cation radius seems to 
be more important than the effect of the A cation. Nev- 
ertheless, density functional calculations^ have shown 
that the defect formation energies are not simple func- 
tions of the cationic radius, and a significant influence 
of the electronic configuration of the A and B cations is 
observed. 

It is now believed that the radiation response cannot 
be described exclusively in terms of the cationic radius 
ratio and the defect formation energies, and that the 
bond-type must be consideredAiSiiiiiii^ It has been sug- 
gested that the resistance of materials to amorphization 
of a complex non-metallic compound is determined by 
a competition between short-range covalent interactions, 
and long-range ionic forces This picture is based on 
experimental evidences for more than 100 different ma- 
terials, with very different chemical and structural com- 
positions. Unfortunately, it is not easy to quantify the 
covalency/ionicity and relate this to radiation resistance, 
specially when the topological freedom of the lattice struc- 
ture can play an important roleit 

The broad range of chemical properties that can be ob- 
tained changing the composition of A- and B-site cations 
in the pyrochlore structure, makes them good testing sys- 
tems to study the relationship between electronic struc- 
ture and resistance to amorphization by radiation dam- 
age. Compounds with strong ionic character, like the 
zirconates (A2Zr20r) are known to be more resistant to 
radiation to clearly covalent stannates (A2Sn207). Ex- 
periments with Gd2(Zr 2; Tii_ :E )207 show that the radi- 
ation resistance increases with increasing Zr-content, a 
clear indication of the effect of the electronic structure. 
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Here we present first principles calculations of two fam- 
ilies of pyrochlores: La 2 B 2 07 and Y2B2O7, with B=Ti, 
Zr, Sn or Hf with very different bonding properties. The 
structural properties and elastic constants are compared 
with available experimental information 1 ^, as well as with 
data obtained with classical atomic models. The struc- 
tural information obtained with these ab initio calcula- 
tions can be used to improve the classical interatomic 
potentials available. 

We will discuss the electronic properties of the different 
compounds and try to quantify the ionicity/covalency of 
the different bonds. Static and dynamical charges are 
frequently used to characterize the nature of bonds in 
molecules and solids. Static charges associated to an iso- 
lated atom are intuitive but ill-defined quantities. Dy- 
namical charges, on the other side, are less intuitive, but 
appear as a more fundamental quantity. These dynami- 
cal charges are very sensitive to the structural properties, 
but also to partial hybridization between occupied and 
unoccupied states^ and hence they can give valuable 
information on the nature of the atomic bonds. 

The paper is organized as follows. In Sec. ^ we will 
describe the structural properties of the pyrochlores stud- 
ied. Section II I II analyses the electronic properties of the 
different bonds, in terms of static and dynamical charges, 
charge density distribution, and decomposition of the 
partial density of states (PDOS). In Sec. II VI a discussion 
of our results is presented, and experimental evidences 
for radiation resistance are examined. Finally, conclu- 
sions are given in Sec. [V] 



II. STRUCTURAL PROPERTIES 

The ideal pyrochlore structure is isometric with space 
group Fd3m, and eight molecules in the unit cell. It can 
be visualized in a variety of ways^S, but is often con- 
sidered as an ordered defect fluorite structure. The A 
and B metal cations occupy the 16d (5,5,5) and 16c 
(0, 0, 0) sites respectively, and the oxygens are in the 48f 
(x, |, |) and 8b (|, |, |) positions. The anion sublattice 
could be completed adding the missing oxygen in the 8a 
site to form the fluorite structure. The system is com- 
pletely described in terms of the lattice size a, and x 
parameter defining the position of the 048f. Geometri- 
cally, the smaller B +4 cation is surrounded by six O48/ 
in a distorted octahedron, and the large A +3 cation is in 
a distorted cubic polyhedron formed by two Osh and six 
C>48f. The 48f oxygen is coordinated to two A and two B 
sites, and the 8b oxygen is inside a tetrahedron formed 
by A +3 cations. 

We have used a variable cell conjugate-gradient mini- 
mization of the energy following the forces and stresses 
to obtain a relaxed minimum-energy structure for each 
of the studied systems. First principles calculations are 
performed with the self-consistent siesta method^ us- 
ing Density Functional Theory (DFT) 22 within the Lo- 
cal Density Approximation (LDA). 23 Norm-conserving 



pseudopotentials 24 in the Kleinman-Bylander formS 5 . 
arc generated with the following atomic valence 
configurations: Ti(3s 2 ,3p 6 ,3d 2 ,4s 2 ), Y(5s 2 ,4p 6 ,4d 1 ), 
Zr(4s 2 ,4p 6 ,4d 2 ,5s 2 ), Sn(5s 2 ,5p 2 ), La(5s 2 ,5p 6 ,5d 1 ,6s 2 ), 
and Hf(6s 1 ,5d 3 ). These, include nonlinear partial-core 
corrections^ and/or semicore states to account for the 
large overlap between core and valence states. We used 
a single-C basis set for the semicore states and double-^ 
plus polarization for the valence states. The charge den- 
sity is projected on a real space uniform grid with an 
equivalent plane-wave cutoff of 280 Ry, to calculate the 
exchange-correlation and Hartree matrix elements in the 
Hamiltonian. Brillouin zone summations are carried out 
with a 3x3x3 Monkhorst-Pack k-point mesh. 



A. Structural parameters 

The relaxed lattice sizes, and x parameters are given 
in Table [I] In general, the calculated LDA lattice pa- 
rameters are underestimated, but the agreement with 
available experimental dataS&i 3 ^! is better than ~ 1%, 
There is also good agreement with other first principles 
calculations >i2*ii However small deviations exist when x 
is compared to the values calculated using classical mod- 
els. These calculations predict an increase in x with the 
B cation radius^. We do not observe this clear trend, 
although the differences appear mainly for the stannates, 
where the covalent character of tin makes the definition 
of ionic radii, at least, questionable. To further validate 
the calculated structures for these pyrochlores, we show 
111 Tab|H the most relevant bond distances and angles. 
The agreement again is remarkable with deviations to 
experimental values not bigger than ~ 2%. 



B. Elastic constants 

The linear elastic constants can be obtained from first 
principles with the derivatives of the stress as a function 
of a properly chosen lattice distortion 6 parametrizing 
the strain^ For a cubic structure, the number of inde- 
pendent elastic constants is reduced to three: en, ci 2 and 
C44. They can be obtained through volume compression, 
tetragonal and trigonal strains: 30 
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where B and C are the bulk modulus and the shear 
modulus. When the applied strain reduces the symme- 
try of the crystal structure adding new internal degrees of 
freedom for the atomic positions, these have to be fully 
relaxed to obtain the elastic constants. In Table [H] we 
show the results for the elastic constants and bulk mod- 
ulus, both for the distorted cells without geometry relax- 
ations (c^- ) , and for those structures resulting in internal 
atomic relaxations. It can be seen that the relaxation 
has a strong effect on the elastic constants, specially for 
c 44 that is softened up to a 40% when the atomic posi- 
tions are optimized. To the best of our knowledge, no 
experimental data is available for the pyrochlores stud- 
ied in this work. It has to be considered that the values 
presented in Table HT1 were obtained at T = OK, and that 
temperature effects generally reduce the values for the 
elastic constants. Consequently, we expect the experi- 
mental values at room temperature to be smaller than 
the values presented here. 

From the elastic constants we can obtain some infor- 
mation on the thermal properties of the pyrochlores. The 



TABLE I: Calculated lattice parameters for A2B2O7 py- 
rochlores, compared to experimental and calculated values. 
Th.A shows results for other first principles simulations. Th.B 
shows results obtained with empirical potentials. 
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ds-o 48f 
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La2Sn207 








This work 


10.744 


0.3341 2.33 


2.61 


2.10 


129.1 


Exp.[2JJ 


10.703 


0.329 2.32 


2.63 


2.07 


131.6 


Th.B[7] 




0.324 












La2Zr2 07 








This work 


10.724 


0.3314 2.32 


2.62 


2.09 


130.5 


Exp. [8] 


10.737 


0.331 








Exp. [7] 


10.805 


0.333 2.34 


2.64 


2.10 


130.6 


Th.A[10l 


10.986 


0.330 








Th.B[7] 




0.326 












La 2 Hf 2 7 








This work 


10.673 


0.3298 2.31 


2.62 


2.07 


131.4 


Exp. [8] 


10.750 


0.331 












La 2 Ti 2 07 








This work 


10.366 


0.3227 2.24 


2.59 


1.98 


135.6 


Th.A[10l 


10.541 


0.323 2.28 


2.63 


2.02 


135.1 






Y 2 Sn 2 7 








This work 


10.348 


0.3437 2.24 


2.44 


2.07 


124.1 


Exp. [27] 


10.372 


0.337 2.25 


2.49 


2.04 


127.6 


Th.A[ll] 


10.329 


0.338 








Th.B[7] 




0.329 












Y 2 Zr 2 7 








This work 


10.335 


0.3417 2.24 


2.45 


2.06 


125.1 


Th.A [11] 


10.463 


0.342 2.26 


2.48 


2.08 


125.0 






Y2H12O7 








This work 


10.300 


0.3399 2.23 


2.46 


2.04 


126.1 






Y 2 Ti 2 7 








This work 


9.974 


0.3321 2.16 


2.43 


1.94 


130.2 


Exp. [13] 


10.100 


0.330 2.18 


2.46 


1.96 


131.3 


Th.A[ll] 


10.049 


0.329 








Th.B[7] 




0.325 









TABLE II: Calculated elastic constants (in GPa) for A2B2O7 
pyrochlores. 
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La 2 Ti 2 7 


336 
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309 


211 


241 


196 


213 


211 


Y 2 Ti 2 7 


388 


160 


318 


236 


381 


153 


241 


229 


La2Zr20 7 


313 


164 


305 


213 


290 


156 


200 


200 


Y 2 Zr 2 7 


394 


170 


312 


245 


351 


163 


198 


225 


La 2 Hf 2 7 


339 


166 


327 


224 


322 


149 


214 


207 


Y 2 Hf 2 7 


403 


191 


343 


262 


379 


167 


221 


238 


La2Sn20 7 


295 


131 


280 


185 


295 


131 


183 


185 


Y 2 Sn 2 7 


354 


141 


289 


212 


348 


140 


202 


209 



longitudinal and transverse acoustic velocities (along the 
[100] direction) are estimated with vi = \fcy\Jp and 
i't = v / C44/ 'p, where p is the mass density of the ma- 
terial^ We can then define the mean sound velocity as: 



We observe that the resulting velocities (shown in Fig^) 
decrease with the atomic mass of the B cation, and also 
with the mass of A. We do not observe a clear trend 
with the cationic radii. These velocities can also be 
used to calculate Debye's characteristic temperature 32 
e D = (hv/k B )(3N /4TrV) 1/3 , for a solid with N Q atoms 
in the volume V, and where h is Plank's constant, fcs is 
Boltzmann constant (see inset in Fig^l. 



III. ELECTRONIC PROPERTIES 

The concept of atomic charge is widely used in chem- 
istry as well as in solid state physics. Unfortunately, it 
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FIG. 1: Calculated mean sound velocity as a function of the 
atomic mass of the B site for the two families of pyrochlores 
studied: Y2B20 7 (empty symbols) and La2B20 7 (filled sym- 
bols). The Debye temperatures O (in K) are shown in the 
inset. 
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is not a well defined concept. Many different definitions 
have been proposed and there is not a general agree- 
ment in how to partition the electronic charge density 
around the atoms and quantify the atomic charge. Be- 
hind this technical problem lies a more fundamental one: 
the quantification of the covalency/ionicity of chemical 
bonds. Here we present the results obtained with several 
definitions of atomic charges, in order to characterize the 
nature of the different bonds in the family of pyrochlores. 

A. Dynamical charges 

The so called dynamical charges (or Born charges) are 
defined in terms of the change of the polarization cre- 
ated when an atom is displaced from its equilibrium po- 
sition. They have been used to characterize the ionic- 
ity/covalency in ferroelectric perovskites^ The tensorial 
definition is: 



where Pp is the (3 component of the macroscopic polar- 
ization induced per unit cell when the atom i is displaced 
in the direction u a , and f2o is the unit cell volume. 

The Born effective charges for the four inequivalent lat- 
tice points are shown in Table ITTT1 The charge neutrality 
sum rule Z* = 0) is frequently used to benchmark the 
accuracy of the calculation, and in our case is fulfilled to 
within 0.04e. The symmetry of the A and B sites allow 
for a decomposition of the Born charge tensors into two 
independent eigenvalues, one for displacements pointing 
in the direction [111] and the others (degenerate) for dis- 
placements in the orthogonal plane. For the A site, the 
highest eigencharge corresponds to the [111] direction, 
due to the presence of the Osb along this axis. On the 
other side, the [111] gives the smallest eigencharge be- 
cause the Oga site is unoccupied in pyrochlores. Small 
deviations from the two-fold degenerate eigencharges are 
observed for some systems, due to the unconstrained ge- 
ometry relaxation. For the O48/ position, the diagonal- 
ization of the symmetric part of the Born tensor gives 
three non-degenerate eigenvectors, along the [Oil], [100] 
and [Oil] directions (with increasing eigenvalues in the 
same order). If we consider the plane defined by the 
B-O-B elements, the first direction corresponds to the 
C2 axis, the second is perpendicular to the plane, and 
the third goes from B to B. Finally, the Ogt, tensor is 
isotropic, and only the diagonal component is given in 
the table. 

The charge for the A site (generally higher than the 
nominal +3e, both for the La and Y pyrochlores stud- 
ied here) is almost independent of the B site compo- 
sition, with variations smaller than 6%. Similar varia- 
tions (~10%) are observed for the Osb- A larger chem- 
ical dependence is observed for the dynamical charges 
associated to the octahedron formed by the B cation and 
the 048f. In this case, there are differences in the Born 



charges for the B cation of up to 40%, and up to 30% for 
the 048f, when the B-site is occupied by Sn instead of Zr 
or Ti. 

The dynamical charges Z* B increase in the order B=Sn, 
Hf, Zr, Ti. The large anomalous charge contributions 
(additional charges to the nominal ionic values of A +3 , 
B +4 and O -2 ) originate from the electronic charge reor- 
ganization induced when an atom is displaced from its 
original position. They are correlated to charge transfers 
between atoms and dynamical changes in the hybridiza- 
tion. Large values of Born charges have been observed 
in Ti-perovskites and attributed to the hybridization be- 
tween the occupied O 2p orbitals and the unoccupied Ti 
3c? orbitalspi^. Hence, it would be naive to conclude that 
the degree of ionicity in pyrochlores increases in the same 
direction (Sn, Hf, Zr, and Ti) as Z* B . 

Figure [21 shows the density of states close to the Fermi 
level for the four lanthanum structures studied (those for 
yttrium- pyrochlore are very similar). The valence band 
has a main contribution from the 2p orbitals of oxygen, 
and a small contribution from the d orbitals of the B 
cation (5s and 5p orbitals for Sn). The conduction bands 




E(eV) 

FIG. 2: (Color online) Electronic density of states (DOS) for 
the La2B2C>7 pyrochlores, with B=Sn, Ti, Hf and Zr (from top 
to bottom panel). The top of the valence band is at eV. The 
dashed, and dot-dashed lines show the projected contribution 
into the C>2p and Lasd atomic orbitals. The thin solid line 
shows the projection into the Sn 5s and 5p orbitals, the Ti 
3d, Zr Ad, and Hf 5d, for each system. 
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TABLE III: Calculated dynamical effective charge tensors for the A ion centered at (|, j, |), B ion at position (0,0,0), 04g/ 
ion at (x, |, |), and Ogi, ion at (|, |, |). The effective charge tensor for any other atom can be obtained applying the symmetry 
operations of the crystal. The eigenvalues of the symmetric part of the tensor are also given. 
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have a sharp peak with B cation d orbital character in 
its lower part, and contributions from the bd orbitals of 
the La atom. La2Sn2C>7 has a very dispersive conduction 
band with tin-5s character. The anomalous contributions 
for the B cations can be interpreted as a result of electron 
current flowing from the occupied 02 P states into the 
empty states of the B-cation in the conduction band when 
the B ion is displaced from its equilibrium position. The 
calculated band gaps are 2.6, 2.8, 4.7 and 4.5 eV for 
Sn, Ti, Hf, and Zr respectively. The LDA band gaps 
are known to be underestimated, and a correction of the 
strong correlation effects is likely to increase the band 
gap and reduce the Born effective charges^. 



B. Static charges 

Many definitions of atomic static charges can be found 
in the literature. They differ in the way the space is par- 
titioned and the charge is associated to each atom. Prob- 
ably the most popular method is the Mulliken population 
analysis^, where the atomic basis functions are used to 
associate the electronic density around each atom. Un- 
fortunately, it is well known that the atomic charges ob- 
tained with this method are basis-dependent. Recently, 
Fonseca-Guerra et. al. showed that the Voronoi Defor- 
mation Density (VDD) method give chemically meaning- 
ful charges with basis set independenc}^. This method 
is based in direct integration of the deformation den- 
sity Sp(r) = p{r) — p a tom(i") over the Voronoi polyhe- 
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portant dependence on the chemical composition appears 
in the charge of the cation B site, with increasing values 
in the order Sn, Ti, Hf and Zr. If the degree of ionicity 
of the pyrochlores is to be associated to the static atomic 
charges this very ordering would be the order of increas- 
ing ionicity (larger atomic charge). Notice that the O48/ 
anions also accumulate more negative charge as we move 
from Sn to Zr, showing that the B-O48/ pairs become 
more ionic. To further validate this assignment, we also 
show the overlap Mulliken population between the oxy- 
gen atoms and the atoms in the A and B sites^ (lower 
panel). These values, though basis-dependent, can be 
used to estimate the degree of covalency (the larger the 
overlap, the larger the covalency). It is shown that the 
overlap population of the oxygen decreases in the order 
Sn>Ti>Hf>Zr. 



C. Charge distribution 

Figure 0] plots the contour of electronic charge re- 
distribution Sp(r) in a plane defined by A-B-O48/, f° r 
Y 2 Ti 2 7 , Y 2 Hf 2 7 , and Y 2 Zr 2 7 . The valence config- 
uration for Ti, Hf and Zr is similar, as well as the lat- 
tice structure. Consequently, the differences in the con- 
tour plots have to be attributed to the difference bonding 
character with oxygen in the three systems. The stan- 
nate is not included here, because of the different valence 
configuration of Sn. The directional character of the Ti— 
O bond is in contrast with the mainly ionic Zr-O. The 
hafhate presents an intermediate directionality. Similar 
contour plots are observed for the lanthanidcs. Notice 
that the charge distribution around the A site is mainly 
spherical showing the high ionicity. The contour plots 
support the hypothesis of increasing covalency from Ti 
to Hf to Zr based pyrochlores. 



FIG. 3: Static charges calculated with VDD (squares) and 
Hirshfeld (circles) methods. Upper panel: the cationic 
charges for the A and B sites (empty and filled symbols). 
Middle panel: oxygen 48f (filled) and 8b charges (empty). 
Lower panel: Overlap Mulliken populations between B and 
O48/ for the two families studied. 



dron around each atom. The charges obtained with 
VDD are similar to the values obtained with the Hir- 
shfeld method^, where the electron density is assigned 
to each atom by weighting with the isolated atom density 

Patom{r)- 

In FigEl we show the calculated static charges with 
VDD and Hirshfeld methods for the cations A and B, 
and for the two anion sites 8b and 48f. We observe that 
the values of Hirshfeld and Voronoi charges are very sim- 
ilar. The composition of the A site affects slightly the 
values for other atoms, specially for the Osb- This is 
because the 8b oxygens are surrounded by just four A 
cations, and no B cation is close enough. The most im- 



IV. DISCUSSION 

Understanding response of materials to radiation is a 
very complex problem. Many factors have been named to 
be relevant in the context of resistance to amorphization 
by irradiation^, including topological freedom, glass- 
forming ability, melting and crystallization temperatures, 
ionicity, bond energy, elasticity, ratio of ionic radii, de- 
fect formation energies, etc. Many criteria to character- 
ize radiation resistance focus in a particular subset of 
properties. In pyrochlores, the defect-formation energies 
have been a popular criterion. It was thought that the 
ionic radius or ionic radius ratio were the most important 
factors in determining the defect energies^ but recent 
first-principles calculations have demonstrated that the 
electronic structure plays an important roleiii The fact 
that all pyrochlores have the same structure, makes them 
an excellent playground to tests the relevance of some of 
these factors, and particularly the role of the bond-type. 

Naguib and Kelly proposed a criterion for resistance 
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FIG. 4: (Color online) Electronic charge redistribution in the A-B-O48/ plane. The lines are taken at intervals of 0.005e in the 
range -0.25c to 0.25e. Solid (red) and dashed (blue) lines correspond to positive and negative values respectively. 



to amorphization by irradiation based on the good em- 
pirical correlation between ionicity and susceptibility to 
amorphization^. Frequently, the ionicity of the material 
is quantified from the difference in Pauling or Phillips 
electro-negativities for different atoms resulting in in- 
consistencies of the Naguib and Kelly criterion. Nev- 
ertheless, the importance of ionicity/covalency in the 
radiation resistance is evident. For the family of py- 
rochlores studied here, the static atomic charges, and 
overlap Mulliken populations indicate that ionicity in- 
creases as the B site is occupied by Sn, Ti, Hf and Zr. 
Remarkably, this is consistent with the increase of resis- 
tance of these materials. The critical temperature (T c ) 
above which complete amorphization of the material can 
not be obtained is frequently used as a measure of the re- 
sistance to amorphization of the material. The lower T c , 
the more resistant the material. For the family of lan- 
thanide pyrochlores La 2 Sn 2 07 (1025 K), La 2 Hf 2 07 (563 
K) and La 2 Zr 2 0y (339 K), the trend agrees with the de- 
crease of covalency, and with the Naguib and Kelly crite- 
rion. Furthermore, the observed^ increase in resistance 
to amorphization with increasing zirconium content in 
Gd 2 Ti 2 _ I Zr x 07 seems to follow the same trend. How- 
ever, Y 2 Sn 2 C>7 with a critical temperature similar to the 
one observed for the hafnate, deviates from this trend 
and would require further study. 

The deformation density contour plot for A 2 Ti 2 07 in- 
dicates that these pyrochlores have an important covalent 
character, as opposed to A 2 Hf 2 07 and A 2 Zr 2 C>7 where a 
more ionic character can be deduced from the deforma- 
tion density. The nearly identical defect formation ener- 
gies calculated for the Ti and Sn members in reference 
could be a consequence of the similar covalent bond. 
We believe that the highly anomalous Born charge cal- 
culated for Ti has an origin on the hybridization between 
the 2p orbitals of oxygen and the 3c? orbitals of Ti, simi- 



larly to the anomalous Born charges observed for example 
in Ti-perovskites, or Ti0 2 , and should not be attributed 
to a strongly ionic character of these pyrochlores. 

The calculated elastic constants do not satisfy the 
Cauchy relation cyijc^ — 1. This condition has to 
be satisfied in crystals with two-body interactions, as 
in a purely ionic system. In particular, for La 2 Sn 2 07, 
La 2 Hf 2 07, Y 2 Ti 2 C>7 and Y 2 Sn 2 C>7 we have C12/C44 < 
0.75. The presence of non-central forces, such as those 
coming from covalent bonding produces deviations from 
Cauchy equality. The relation is more strongly broken for 
the unrelaxed elastic constants (c^/cjj.j, and is partially 
recovered under relaxation. The surprising recovery ob- 
served for La 2 Ti 2 C>7 could be related to the instability of 
the pyrochlore structure for this system. 



V. CONCLUSIONS 

In this work we have used density functional theory 
calculations to analyze the structural and chemical prop- 
erties of a family of pyrochlores. The lattice and struc- 
tural parameters (a and x) were obtained with a conju- 
gate gradient minimization of the forces and strains. Ap- 
plying a particular set of lattice distortions, the elastic 
constants were also obtained. The structural properties 
obtained from first principles can be used to parametrize 
new interatomic potentials that could be used in classi- 
cal Molecular Dynamics simulations for materials design. 
The lack of an accurate description of the chemical bond 
in these classical simulations was criticized recently in the 
context of resistance to amorphization by irradiation. 

Although experimental evidence shows that the elec- 
tronic structure has to play some role in radiation resis- 
tance, it is still not clear what is the trend. The similar 
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lattice structure for the ideal systems (not all of them 
stable in the pyrochlore structure) studied here, makes 
them excellent systems to elucidate the relevance of the 
chemical properties. If the description of Trachenko and 
collaborators is correct, the competition between short 
(covalent) and long (ionic) interactions is strongly related 
to the resistance to amorphizationi5ii£. 

In this work we present a study of some indicators 
of covalency/ionicity in pyrochlores. Although some of 
these parameters (static charges, deformation density) 
seem to show evidence that the stannates and titanates 
pyrochlores have a more covalent character than the 
highly ionic hafnates and zirconates, (which are known 
to be more resistant to amorphization), others, like the 
dynamical charges, do not show the same trend. Fur- 



ther work is required to understand these differences. In 
particular, a systematic experimental study of T c is desir- 
able, as most of the available data is recorded in different 
conditions (bombarding ions, energy, etc.). 



Acknowledgments 

This work was supported by BNFL and UK's NERC 
through the eminerals consortium. The authors are 
grateful to G. Lumpkin, K. Whitlc, K. Trachenko, I. 
Farnan, and S. Rfos for fruitful discussions on radiation 
resistance of pyrochlores. EA thanks the Donostia Inter- 
national Physics Center (DIPC) for its hospitality. 



1 D. P. Carm, C. A. Randall, and T. R. Shrout, Solid State 
Commun. 100, 529 (1996); T. R. Shrout and S. L. Swartz, 
Mater. Res. Bull. 18, 663 (1983). 

2 R. Vassen, X. Cao, F. Tietz, D. Basu, and D. Stover, J. 
Am. Ceram. Soc. 83, 2023 (2000); G. Suresh, G. Seeni- 
vasan, M. V. Krishnaiah, and P. S. Murti, J. Nucl. Mater. 
249, 259 (1997) 

3 B. J. Wuensch, K. W. Eberman, C. Heremans, E. M. Ku, 
P. Omnerud, E. M. E. Yeo, S. M. Haile, J. K. Stalick, and 
J. D. Jorgensen, Solid State Ionics. 129, 111 (2000). 

4 R. C. Ewing, W. J. Weber, and J. Lian, J. Appl. Phys. 95, 
5949 (2004). 

5 K. E. Sickafus, L. Minervini, R. W. Grimes, J. A. Valdez, 
M. Ishimaru, F. Li, K. J. McClellan, T. Hartmann, Science 
289, 748 (2000). 

6 P. K. Schelling, S. R. Phillpot, and R. W. Grimes, Phil. 
Mag. Lett. 84, 127 (2004). 

7 Y. Tabira, R. L. Withers, L. Minervini, and R. W. Grimes, 
J. Solid State Chem. 153, 16 (2000). 

8 G. R. Lumpkin, K. R. Whittle, S. Rios, K. L. Smith, and 
N. J. Zaluzec, J. Phys.: Condens. Matter. 16, 8557 (2004). 

9 M. Pirzada, R. W. Grimes, L. Minervini, J. F. Maguire, 
and K. E. Sickafus, Solid State Ionics 140, 201 (2001). 

10 A. Chartier, C. Meis, W. J. Weber, and L. R. Corrales, 
Phys. Rev. B 65, 134116 (2002). 

11 Wendy R. Panero, Lars Stixrude, and Rodney C. Ewing, 
Phys. Rev. B 70, 054110 (2004). 

12 J. Lian, R. C. Ewing, L. M. Wang, and K. B. Helean, J. 
Mater. Res. 19, 1575 (2004). 

13 J. Lian, J. Chen, L. M. Wang, R. C. Ewing, J. M. Farmer, 
L. A. Boatner, and K. B. Helean, Phys. Rev. B 68, 134107 
(2004). 

14 K. Trachenko, M. Pruneda, E. Artacho, M. T. Dove, Phys. 
Rev. B 70, 134112 (2004). 

15 K. Trachenko, J. Phys.: Cond. Matt. 16, R1491 (2004). 

16 K. Trachenko, J. M. Pruneda, E. Artacho, and M. T. Dove, 
Phys. Rev. B 71, 184104 (2005). 

17 L. W. Hobbs, Nucl. Instrum. Methods B 91, 30 (1994). 

18 Not all the compounds in the families studied are stable in 
the pyrochlore structure. For example La2Ti2 07 has the 
layered perovskite structure. 

19 Ph. Ghosez, J. -P. Michenaud, and X. Gonze, Phys. Rev. 
B 58, 6224 (1998). 

20 M. A. Subramanian, G. Aravamudan, G. V. S. Rao, Prog. 



Solid. State Chem. 15, 55 (1983). 

21 P. Ordejon, E. Artacho, and J. M. Soler, Phys. Rev. B 53, 
R10441 (1996); J. M. Soler, E. Artacho, J. D. Gale, A. 
Garcia, J. Junquera, P. Ordejon and D. Sanchez- Portal, J. 
Phys.: Condens. Matter 14, 2745 (2002). 

22 P. Hohenberg and W. Kohn, Phys. Rev. 136, B864 (1964); 
W. Kohn and L. J. Sham, Phys. Rev. 140, A1133 (1965). 

23 D. M. Ceperley and B. J. Alder, Phys. Rev. Lett. 45, 566 
(1980). 

24 N. Troullier and J. L. Martins, Phys. Rev. B 43, 1993 
(1991). 

25 L. Kleinman and D.M. Bylander, Phys. Rev. Lett. 48, 1425 

(1982) . 

26 S. G. Louie, S. Froyen, M. L. Cohen, Phys. Rev. B 26, 
1738 (1982). 

27 B. J. Kennedy, B. A. Hunter, and C. J. Howard, J. Solid 
State Chem. 130, 58 (1997). 

28 L. Minervini, R. W. Grimes, Y. Tabira, R. L. Withers, K. 
E. Sickafus, Philos. Mag. 82, 123 (2002). 

29 O. H. Nielsen and R. M. Martin, Phys. Rev. Lett. 50, 697 

(1983) . 

30 R. Stadler, W. Wolf, R. Podloucky, G. Kresse, J. 
Furthmiiller, and J. Hafner, Phys. Rev. B 54, 1729 (1996). 

31 H. Ibach and H. Liith Solid-State Physics. An introduction 
to Principles of Materials Science ( Springer- Verlag, Berlin 
(1996). 

32 R. A. Robie and J. L. Edwards, J. Appl. Phys. 37, 2659 
(1966). 

33 M. Posternak, R. Resta, and A. Baldereschi, Phys. Rev. B 
50, R8911 (1994). 

34 A. Filippetti and N. A. Spaldin, Phys. Rev. B 68, 045111 
(2003). 

35 R. S. Mulliken, J. Chem. Phys. 23, 1833 (1955). 

36 C. Fonseca Guerra, J.-W. Handgraaf, E. J. Baerends, F. 
M. Bickelhaupt, J. Comput. Chem. 25, 189 (2003). 

37 F. L. Hirshfeld, Theoret. Chim. Acta 44 129 (1977). 

38 G. R. Lumpkin, S. Rfos, K. R. Whittle, J. M. Pruneda, K. 
Trachenko, K. L. Smith, and N. J. Zaluzec, to be published. 

39 H. N. Naguib and R. Kelly, Radiat. Eff. 25, 1 (1975). 

40 S. X. Wang, B. D. Begg, L. M. Wang, R. C. Ewing, W. 
J. Weber, and K. V. G. Kutty, J. Mater. Res. 14, 4470 
(1999). 



